Paired composite fermion wavefunctions 
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We construct a family of BCS paired composite fermion wavefunctions that generalize, but remain 
in the same topological phase as, the Moore-Read Pfafhan state for the half-filled Landau level. It 
is shown that for a wide range of experimentally relevant inter-electron interactions the groundstate 
can be very accurately represented in this form. 



The nature of the fractional quantum Hall effect at 
v = | has been the subject of continued interest since its 
discovery roughly two decades ago [l[. The Moore-Read 
Pfaffian wavefunction Q, describing p-wave pairing of 
composite fermions is currently the best candidate 
wavefunction for this state [H, 0]- Due to the remark- 
able property of having quasi-particles with non-abelian 
statistics, this state has recently attracted interest in the 
context of fault-tolerant topological quantum computa- 
tion [f|. Though the Moore-Read state is a well estab- 



lished candidate for the groundstate at v 



its overlap 



with the exact groundstate in simulations of small sys- 
tems is rather low in comparison to other known trial 
states at different filling factors 0, 0. This is particu- 
larly discomforting as no explicit construction for pertur- 
bations around the Moore- Read state has been previously 
available, and the Moore-Read state has been described 
as a somewhat unique choice [i| . Furthermore several re- 
cent studies have altogether challenged the view of v = § 
as being the Moore- Read state [7J, |8| . 

In this Letter, we introduce a general representation 
of paired composite fermion (CF) states, merging the 
concept of BCS Hall states M with the explicit construc- 
tion of CF wavefunctions [lOj, [llj . The Moore-Read state 
can be cast very accurately in this form, which reveals 
its connection to the pairing of composite fermions on 
top of a Fermi-sea, and shows how our general paired 
CF-BCS wavefunctions are adiabatically connected to 
Moore-Read. We also compare our trial states to the 
exact groundstates of the Coulomb Hamiltonian He for 
electrons in the 1 st excited Landau level (ILL), plus an 
arbitrary additional pscudopotcntial 8V\ interaction. For 
a very broad range of 8V\ we find very high overlap of 
our trial wavefunctions with the exact groundstate, thus 
showing the extent of the Moore-Read phase. 

For our description of the physics at v = |, we shall 
assume the lowest Landau level (LLL) to be entirely filled 
and inert, such that the relevant degrees of freedom cor- 
respond to a half filled ILL. We assume the spin degree 
of freedom of these electrons to be frozen (although the 
experimental situation is less clear [13]). The ILL is rep- 
resented by wavefunctions in the LLL using appropriately 
modified pseudopotential coefficients [l3| . 

The aim of our construction is to "composite fermion- 
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simple BCS state. In second quantized nota- 



the general form of the BCS groundstate is [14 1 

bcs) = Ilk (l + 9k e itp 4: c -k) 1°) i written in an un- 
normalized manner here. This wavefunction can be pro- 
jected to a fixed number of particles by integration over 
/ d(p exp(-iNip) such that we retain exactly TV pair cre- 
ation operators. The (inverse) Fourier transform into real 
space then yields |14| 



*BCs(n,...,rAr)=?Pf[ff(ri-r i )] 



(1) 



where the Pfaffian ^Pf is an antisymmetrized sum over 
all possible pairings SPf(sij) = ./%i2034 • • -9n-i,n] = 
±aV| det gij\ with A the antisymmetrization operator. In 
Eq. [U g is constrained to be an antisymmetric function, 
given in terms of its Fourier components by 



X> e<k ' (ri " r,) E^Sk^fr^-xfe). (2) 

k k 



For the last equivalence, we have identified the exponen- 
tial factor as the product of two basis functions <^k(r) = 
exp(ik ■ r) of free electrons on the plane. This product of 
free wavefunctions form is naturally generalized to spher- 
ical geometry below. 

In order to obtain a LLL wavefunction at filling fac- 
tor v = i , we follow Jain's approach [lo| of multiplying 
a bare electron wavefunction with Jastrow factors, and 
projecting the result to the LLL, yielding [15 1 



, z N ) 



where "P L ll is the LLL projection operator, and Z\ is 
the complex representation of r.;. The special case g = 
l/(zi — Zj) reproduces the Moore-Read wavefunction (and 
the projection then becomes trivial). 

In order to render the projection Plll numerically 
tractable in general, we bring single particle Jastrow fac- 
tors Ji = Y[k=£i( z i ~ z k) inside the Pfaffian on every line i 
and every column j of the matrix g, ■ , and project each of 



the matrix elements individually 11]. As demonstrated 
in Ref. [llj], such slight changes in the projection prescrip- 
tion do not alter the accuracy of the composite fermion- 
ization procedure. It is thus expected that projecting ma- 
trix elements [</(r,- — r ,•)./,•</,•] individually is very similar 
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to a projection of the full wavefunction. For further sim- 
plification, one can decompose g analogous to Eq. [2l and 
apply the projection separately to each of the orbitals 4>k 
as suggested in Ref. . Again this slight change in pro- 
jection prescription is not expected to damage our wave- 
function. Denoting 4>w(zi) = Jf 1 Vi jL i J [(j)] l (zi)Ji], Jastrow 
factors may be factored again outside the Pfaffian, and we 
obtain the final expression for general composite fermion- 
ized BCS (CF-BCS) states 



J^5k 4>k(Zi) 4>-k(Zj) 
k 



(3) 



In the remainder of this study we will focus on finite 
size svstems with N electrons on the spherical geome- 
try [13]. In order for Eq. [3] to represent the Moore- 
Read phase, we must work at a flux of N$ = 2N — 3. 
The orbitals </>k thus correspond to composite fermions 
in one quantum of negative effective flux (ljj . ie., the 
(very small) effective magnetic field experienced by CF's 
is directed opposite to the external magnetic field. The 
relevant CF orbitals (</>k) are given by the projected 

monopole harmonics Yn~m 2 studied in Ref. 0. To 
assure @ that the angular momentum of a pair is 
I = — 1 (the negative p-wave pairing of the Moore- 
Read phase) we must choose g^ — > (— l) q+m g n and we 
are left with only one variational parameter g n per CF 
shell. Thus the term in the brackets of Eq. [3J becomes 
En, m (^r- 1/2 9nYr m ~Hz l )Y:=-J(z 3 ). The sum over 
n goes from n — Qton — N~2 since orbitals with 
n > N — 1 are projected to zero. Thus, up to an overall 
normalization, there are N — 2 variational parameters. 

It is also possible to study other pairing symmetries 
within our approach. These would yield states at differ- 
ent values of the flux N$. Here, we focus on negative 
p-wave pairing, which appears most consistent with pre- 
vious numerical data 0, HJ ■ 

The variational character of the wavefunctions we 
study (Eq. [3]) implies that we need to optimize over the 
set of parameters g = (go, <?i,S2, ■ ■ ■ ,9n-2) to obtain a 
good trial wavefunction. The definition of a "good" wave- 
function is somewhat arbitrary, and one may attempt to 
optimize various measures of its accuracy, e.g. the energy 
of the wavefunction, the overlap with the exact ground- 
state, or the error in the pair correlation function com- 
pared to the exact groundstate. The chosen measures of 
accuracy are evaluated by Monte-Carlo, and the varia- 
tional parameters are then optimized. 

It is instructive to verify that the Moore-Read state can 
be reproduced as a CF-BCS state (Eq. [3]). Numerically 
we find that for a suitable set of variational parameters, 
g, we are able to achieve overlaps in excess of 0.99 with 
the Moore- Read state for systems with up to 20 electrons. 
While this may seem a rather complicated reformulation 
of the Moore-Read state, we can now perturb the wave- 
function with our variational parameters. 



Fig. [T] shows overlaps between trial states and the cor- 
responding exact groundstates for different interactions 
parametrized by the pseudopotential 5Vi (with SVi = 
being the pure Coulomb interaction in the ILL, and very 
large 5V\ roughly corresponding to the interaction of the 
LLL). We also show overlaps of the exact groundstate 
with the Moore-Read state and the CF liquid (CFL) 
(here defined to be Eq. [3J with the occupation coeffi- 
cients g n being unity below the chemical potential and 
zero above). Results are shown for TV = 12, 14, 16. The 
dimensions of the L — Hilbert space are respectively 
d = 52,291,2077. Although we could in principle opti- 
mize over as many as N — 2 variational parameters, in 
practice we find very good wavefunctions using at most 
the first 7 parameters. We remind the reader that the 
variational parameters are included (like the it's and v's 
of BCS theory) to optimize the shape of the pairing wave- 
function. Note that the number of parameters used is far 
less than the dimension of the L = Hilbert space, so the 
good agreement with exact diagonalization is significant. 
Further we emphasize that the trial states (Eq. [3j have 
very high overlap with the exact groundstate for a wide 
range of values of 8 V\ . 

For a particular value of 8V\ ~ 0.04 the Moore-Read 
state is also a very good trial state. However even at 
this value of 6Vi , our trial states yield an improved rep- 
resentation of the exact groundstate. For larger 8V\ the 
CF-BCS states become even more accurate and can be 
continuously deformed into the composite fermion liquid 
at large 8V\ . As shown in Fig. Q] when 8V\ gets close to 
zero (or negative, not shown) the overlap drops. This 
behavior could be expected considering prior work Q 
showing a nearby phase transition, as we discuss below. 

As mentioned above, overlaps are not the only possible 
measure of the accuracy of a wavefunction. In Fig. [2] we 
show how the pair correlation functions h{&) of different 
trial states compare against those of the exact ground- 
state. Here 6 is the angle between two particles on the 
sphere, and we show the mean square error in the correla- 
tion function Sh 2 = J d(cos6>) \h(0) — h cxact (9)\ 2 . We note 
that if 8h 2 — for a trial state, then it is identical to the 
exact groundstate (this can be seen from the variational 
principle, noting that h{6) fully determines the energy 
for a two-body interaction) . Fig. [5] once more illustrates 
that our trial wavefunctions are extremely accurate - far 
more so than either the Moore-Read or CFL trial states. 
Again, we find that near 8V\ = our trial state fails to 
match the exact pair correlation function to some extent. 

In Fig. [21 the trial wavefunctions have been reopti- 
mized with respect to 8h 2 . The overlaps of these new 
wavefunctions with the exact groundstate, would gener- 
ally be found to be slightly lower than those in Fig. [IJ but 
still remain very high (except in the vicinity of 8V\ = 0). 

Since, as mentioned above, we can continuously deform 
our wavefunctions to have over 0.99 overlaps with the 
Moore-Read state, we conclude that our CF-BCS wave- 
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FIG. 1: (color online) Overlaps of trial states with exact 
groundstate as a function of the interaction parameter SVi 
for N — 12, 14 and 16 electrons. Here, SVi = corresponds 
to the pure Coulomb interaction in the ILL. The optimized 
composite-fermionized BCS wavefunctions (Eq. [3] black cir- 
cles) have very high overlap except close to SVi = 0, where 
the system is thought to be close to a phase transition. The 
Moore-Read wavefunction (blue squares) is also good near 
8Vi = 0.04, but falls off substantially at other values. The 
composite fermion liquid wavefunctions (orange diamonds) 
are accurate at very high SVi only. Error bars indicate sta- 
tistical errors where a Monte-Carlo algorithm was employed 
for the evaluation of the overlaps. The high accuracy of the 
BCS wavefunctions over a broad range of interactions shows 
the large extent of the weak-pairing phase. 



functions are generally in the same so-called weak-pairing 
phase as the Moore-Read state. To further emphasize 
this point, we note that a wavefunction in a weak-pairing 
phase should have the property [9( that g(r) ~ 1/z at 
large distances r. While this is not obvious from the 
form of Eq. [3] (particularly considering the complexity 
of the projected wavefunctions <fi) we can nonetheless es- 
tablish it is true in several ways. Firstly, we have tried 
making the 1/z tail of the pair correlation function ex- 
plicit, writing gfa —tj) = a/ (zi — Zj) + /(r^ — r 3 ) before 
projection, decomposing only the function / into orbitals 
as in Eq. [5] and projecting these orbitals. We have found 
that this procedure leads to equivalent results. Secondly, 
the property g(r) ~ 1/z at large r implies that the k — » 
orbitals are occupied with probability approaching unity 
Q (which would not be true of a strong pairing phase) . 
It is easy for us to establish numerically that the lowest 
orbitals (n = 0) are indeed fully occupied unity by test- 
ing that increasing the value of the variational parameter 
go does not change the wavefunction. 

It is also worth checking that the exact groundstate is 
indeed adiabatically connected to the Moore-Read state. 
To this end, we analyze the evolution of the energy 
gap for a family of Hamiltonians that interpolate be- 
tween the three-body contact interactions V3 , which yield 
the Moore-Read state as its exact groundstate, and a 
two-body interaction Hamiltonian TL' C corresponding to 
SVi = 0.04. In particular, for any of the interactions 



FIG. 2: (color online) Squared error 5h 2 in the pair corre- 
lation function of various trial wavefunctions compared with 
exact groundstates as a function of the interaction parame- 
ter SVi for N — 12, 14 and 16 electrons. Symbols are as in 
Fig. HI Again we find that composite-fermionized BCS wave- 
functions are far more accurate than either Moore-Read or 
the CF liquid. Inset: Pair correlation functions h(8) of the 
Moore-Read state and our trial wavefunction along with the 
exact groundstate at N = 14 and SVi = 0.02. Our trial state 
is essentially indistinguishable from the exact groundstate, 
whereas the Moore-Read state is slightly different. 

Ti(x) = xV 3 + (1 — x)H' c , we find no indication that in 
the thermodynamic limit the energy gap closes (data not 
shown). We conclude that the exact groundstate of H' c 
is adiabatically connected to the Moore-Read state, con- 
firming that the exact groundstate of SV\ = 0.04 is in the 
weak pairing phase [t| . 

The main result of this work is the construction of a 
family of accurate wavefunctions in the same topologi- 
cal phase as the Moore-Read wavefunction. This can be 
thought of as the composite fermionization of a weakly 
paired BCS wavefunction. We find that over a broad 
range of interactions these wavefunctions are very accu- 
rate - far more so than the Moore-Read wavefunction 
itself, which should be thought of only as an example of 
a wavefunction in a broad phase of matter. Indeed, the 
Moore-Read state may be approximated extremely pre- 
cisely by the form we propose, and when doing so, the 
result does not particularly stand out from other possible 
CF-BCS states. Although from a topological standpoint, 
it is sufficient to identify the phase of matter, from a prac- 
tical standpoint, it is still valuable to have explicit forms 
of wavefunctions [Icl [ll[, as this is important for per- 
forming detailed calculations of excitation spectra and 
other physical properties. Although we have currently 
only analyzed the composite fermionization of ground- 
state BCS wavefunctions, it is natural to consider a sim- 
ilar procedure for excited states, which we will consider 
in future work. 

Let us now discuss how our work reflects on and re- 
lates to previously results. Prior work on the torus Q 
found a first order phase transition from a charge density 
wave (CDW) state to a phase presumed to be the Moore- 
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Read phase at roughly the Coulomb point 8V\ = 0. This 
phase had the required degeneracy of a weak pairing 
(Moore-Read) phase, but had relatively low overlaps with 
the Moore-Read wavefunction itself. Particle-hole sym- 
metrizing the Moore-Read wavefunction increased the 
overlap to 97% for N = 10 at one particular value of 5V\ 
but remained somewhat lower at other values. We note 
that the Moore-Read phase and its particle-hole conju- 
gate are distinct phases [13] and the effect of symmetriza- 
tion is unclear and remains a topic of current interest. In 
our work on the sphere, there is no possible mixing of 
states with their conjugates, although we cannot deter- 
mine whether a state or its conjugate would occur in an 
experimental system. On the sphere, it was previously 
known [3] that the overlap of the exact groundstate with 
the Moore- Read state has a peak at 6V\ « 0.04, and also 
drops strongly near SVy — 0- However, on the sphere, it 
was hard to distinguish the thermodynamic phase since 
there is no groundstate degeneracy to use as a guide. Our 
work, in contrast, studies only trial wavefunctions in the 
Moore-Read phase. 

In contrast to all prior work, our trial wavefunctions 
have high overlaps over a very broad range of 5Vi, con- 
firming that the weak pairing phase is robust to large 
changes in the interaction. Our wavefunctions make a 
smooth transition between the Moore-Read phase and 
the CF liquid at large 5V\. It is difficult to distinguish 
numerically if the groundstate of the LLL still has some 
amount of pairing. To determine if at large SVu the puta- 
tive CFL still pairs (as previously suggested [a]), a more 
careful study of the groundstate for the LLL interactions 
would be required. We note in passing that the wide 
region of intermediate values of 8V\ (between where the 
Moore- Read wavefunction is accurate and where the CFL 
becomes accurate) which we describe extremely well with 
our wavefunctions, could be hard to access with typical 
2DEG samples but could likely be realized using hole- 
doped samples [!§] or graphene (l9j |. 

At interactions SVi < 0.04 our wavefunctions have 
substantially better performance than the Moore-Read 
wavefunction. However, near 8V\ ks our wavefunctions 
do not perform as well as one might hope. This is not 
surprising considering that the groundstate of 8V\ = 
on the torus is a CDW state [5[. However, experiments, 
which see a quantum Hall plateau, do not correspond 
to the pure Coulomb interaction (8Vi — 0) due to finite 
well width effects and Landau-level mixing. It has also 
been noted [B[ that a more realistic interaction puts the 
physical system just slightly on the quantum Hall side of 
the transition. Indeed, it is known experimentally (2pj 
that modifying the electron interaction slightly by tilting 
the field pushes the system from a quantum Hall state 
into a CDW state. Being that pure Coulomb is thought 
to be on the other side of this phase transition, the fact 
that our wavefunctions remain so good is perhaps sur- 
prising. However, one might argue that since the CDW 



is frustrated by the geometry of the sphere we can still 
match the groundstate reasonably well with a sufficiently 
perturbed weak pairing wavefunction, which remains adi- 
abatically connected to the Moore-Read state. 

In Ref. it was suggested that the gapped state near 
the Coulomb point is best constructed within a com- 
posite fermions basis without appeal to the Moore-Read 
wavefunction. Our wavefunction is indeed constructed 
in terms of composite fermions, retains relatively high 
similarity with the exact groundstate, and also remains 
adiabatically connected to Moore-Read. We believe this 
result should put to rest concerns that the gapped phase 
near the Coulomb point is not in the topological phase of 
the Moore- Read state or its particle-hole conjugate fl7j ]. 
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